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We consider various random models (directed polymer, random ferromagnets, spin-glasses) in 
their disorder-dominated phases, where the free-energy cost F{L) of an excitation of length L 
presents fluctuations that grow as a power-law AF{L) ~ with the 'droplet' exponent 9. Within 
the droplet theory, the energy and entropy of such excitations present fluctuations that grow as 
AE{L) ~ AS{L) ~ jj-nl"^ where da is the dimension of the surface of the excitation. These systems 
usually present a positive 'chaos' exponent = ds/2— 9 > 0, meaning that the free-energy fluctuation 
of order is a near-cancellation of much bigger energy and entropy fluctuations of order V^'l"^ . 
Within the standard droplet theory, the dynamics is characterized by a barrier exponent V satisfying 
the bounds 9 d — \. In this paper, we argue that a natural value for this barrier exponent is 

i/i = ds/2: (i) for the directed polymer where ds = 1, this corresponds to ijj = 1/2 in all dimensions; 
(ii) for disordered ferromagnets where ds = d — 1, this corresponds to ^ = [d — l)/2; (iii) for 
spin-glasses where interfaces have a non-trivial dimension ds known numerically, our conjecture 
ip = ds/2 gives numerical predictions in d = 2 and d = 3. We compare these values with the 
available numerical results for each case, in particular with the measure t/j ~ 0.49 of Kolton-Rosso- 
Giamarchi, Phys. Rev. Lett. 95, 180604 (2005) for the non-equilibrium dynamics of a directed 
elastic string. 

I. INTRODUCTION 

The non-equilibrium dynamics of extended objects such as polymers or interfaces in random media has remained 
very controversial over the years. The first question is whether the dynamics is algebraic or logarithmic, or equivalently 
whether there exists a positive barrier exponent ■0 as we now recall. 

A. Logarithmic dynamics with some barrier exponent t/j > 

The activated nature of the non-equilibrium dynamics is natural within the droplet scaling theory proposed both 
for spin-glasses P, [1, [1] and for the directed polymer in a random medium [j (see Section [III below for more details). 
Barriers are then expected to grow as a power-law at large scale 

B[L) - L'l' (1) 

where the exponent ip > Q i& constant in the whole low-temperature phase T < Tc and is a property of the large-scale 
zero-temperature fixed point. The typical time ttyp{L) associated to scale L grows as an exponential 

\nttyp{L) ^ B{L) ^ L"^ (2) 

As a consequence, the non-equilibrium dynamics starting at time t = is expected to involve only logarithmic 
functions of time via the characteristic length-scale L{t) associated to time t 

L{t) (Ini)i (3) 

A sirnple one-particle one-dimensional disordered model where this type of activated dynamics occurs is the Sinai 
model [5| where one particle diffuses in a random Brownian potential of exponent ip = 1/2 : the diffusion is then 
logarithmic with L(t) ~ (Int)^. One-time and two-time properties of the non-equilibrium dynamics can be computed 
at large time via a strong-disorder renormalization procedure that yields asymptotic exact results at large times @ . 
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B. Algebraic dynamics with some dynamical exponent z 

In the alternative scenario of 'algebraic' dynamics, barriers grow at most logarithmically with L, i.e. the exponent 
■0 of Eq. [T] vanishes V' = 0. Time scales and length scales are then related by some dynamical exponent z 

ttyp{L) ^ (4) 
Lit) r^ti (5) 

instead of the logarithmic relations of Eqs [2] and [3l As a consequence, the aging properties of the non-equilibrium 
dynamics involve ratios of times (instead of ratios of logarithms of the times ). 

A simple one-particle one-dimensional disordered model where this type of algebraic aging occurs is the Bouchaud 
trap model Q (see Q for the mean-field version). 



C. Debate on the phase-space structure at large scales 



This question on the value of the barrier exponent ^ amounts more generally to understand the phase space 
structure of polymers or interfaces in random media. Within the droplet theory [H, [3, Q (^^^ Section , static 
properties are governed by low-energy excitations characterized by the 'droplet' exponent 9 and by the dimensionality 
ds of the surface of the excitations, whereas the dynamics involves a priori another exponent ip satisfying the bounds 
9 < t/j < d — 1 0, Whereas the exponents 9 and ds for the statics are known either exactly or numerically in 

various models, the numerical measure of the barrier exponent ^ turns out to be much more difficult from the point 
of view of computational complexity ^9j. As a consequence, the value of "0 has remain very uncertain even numerically 
in many random models. The aim of this paper is to explain that the value -0 = ds/2 for the barrier exponent is very 
natural for disordered models that have a positive chaos exponent — ds/2 — 9 > 0, and to compare with available 
numerical data. 



D. Organization of the paper 



In Section [TTl we recall the essential properties of the droplet scaling theory and present our general arguments 
for the conjecture ip — ds/2. We then discuss this conjecture with a comparison to existing numerical results for 
the following models : directed polymers in random media (Section IIIip . disordered ferromagnets (Section livp . and 
spin-glasses (Section [V]). Section |VT] contains our conclusions. 



II. ARGUMENTS IN FAVOR OF THE VALUE ip = d,/2 FOR THE BARRIER EXPONENT 



In this section, before explaining the conjecture -0 — ds/2 for the barrier exponent, we need to recall the main 
statements of the droplet scaling theory proposed both for spin-glasses [H, [1| and for the directed polymer in a 
random medium M. 



A. Reminder on equilibrium properties within the droplet theory 

1. Statistics of excitations above the ground state 

At very low temperature T — > 0, all observables are governed by the statistics of low energy excitations above the 
ground state. An excitation of large length I costs a random energy 

E,,cil,T = 0)^fu (6) 

where 9 is the so-called 'droplet' exponent [J, 0, 0, 0], and where u is a positive random variable distributed with 
some law Qo{u) having some finite density at the origin Qo{u = 0) > 0. A low-temperature disorder-dominated phase 
exists whenever the droplet exponent 9 is positive 9 > 0. 

From Eq. [6l the probability distribution of large excitations I ^ 1 reads within the droplet theory 

dlpeUl) - ~ fe-P''- (7) 
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where the factor dl/l comes from the notion of independent excitations In particular, its average over the disorder 
foUows the power-law 

/• + 00 11 11 

dl^:jr)^j^ duQo{u)je-f''"- = TQ{0)^ (8) 

Since correlation functions at large distance are directly related to the probability of large excitations, the low 
temperature phase is very non-trivial from the point of view of correlations lengths : the typical exponential decay of 
Eq[7] indicates a finite typical correlation length ^typiT), whereas the averaged power-law behavior of Eq. [5] means that 
the averaged correlation length S,aviT) is actually infinite in the whole low temperature phase ^au(0 < T < Tc) = oo. 



2. Low temperature phase governed by a zero-temperature fixed point 

According to the droplet theory, the whole low temperature phase < T < Tc is governed by a zero-temperature 
fixed point. However, many subtleties arise because the temperature is actually 'dangerously irrelevant'. The main 
conclusions of the droplet analysis 0, Q can be summarized as follows. The scaling of Eq [5] governs the free energy 
cost of an excitation of length I, provided one introduces the typical correlation length £,typ{T) to rescale the length I 



F,Ul, 0<T<T,)^ WTO ) 

Here as before, u denotes a positive random variable distributed with some law Q(u) having some finite density at 
the origin Q{u = 0) > 0. But this droplet free energy turns out to be a near cancellation of much larger energy and 
entropy contributions that scale for large I as 0, 0| 

E,,,{l,0<T <T,) ^a{T)&w + ei{T)f (10) 

where ds represents the dimension of the surface of the excitation. The random variable w of order 0(1) and of zero 
mean is expected to be Gaussian distributed. The argument is that the energy and entropy are dominated by small 
scale contributions of random sign 0, whereas the free energy is optimized on the coarse-grained scale £,typ{T). A 
very important consequence of the difference in scaling of the free-energy fluctuation of Eq[9]and of the energy-entropy 
fluctuations of Eq. [10] is the presence of disorder and temperature chaos in the whole low-temperature phase with the 
so-called chaos exponent 0, Q, [l3| 

C=y-e>0 (11) 

Note that numerically, temperature chaos is usually harder to observe than disorder chaos (see [Til [T3 | and references 
therein). 

For numerical simulations, it is important to stress that the term of order /''"/^ in Eq[TU]is the leading term for 
large I, but that there exists a sub- leading term of order to recover the free-energy fluctuations of Eg El And since 
the amplitude criT) of the leading term vanishes in the limit of T = (as cr{T) oc T x T^/^ — T'^l'^ [ll|), whereas 
the amplitude ei(T) of the subleading term in the energy remains finite e\(T = 0) > 0(Eq. (5]), one needs to simulate 
sufhciently large excitations to reach the size / where the leading term of Eq. [10] becomes much bigger than the 
subleading term : 

a{T)l^ ^ ei{T)f (12) 



B. Reminder on the non-equilibrium dynamical properties within the droplet theory 

Within the standard droplet theory P, @, the non-equilibrium dynamics is governed by large-scale barriers 
B(V) ^ where the barrier exponent satisfies the bounds 

e<^<d-l (13) 

The lower bound comes from the fact that the barrier B{L) to create a droplet excitation of size L cannot be less than 
the free-energy cost of the droplet (Eq[5]). The upper bound comes from the expectation that the barrier cannot have 
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a greater exponent than the barrier L needed to create a non-optimized excitation of surface L . In particular, 
whenever there exists a low-temperature disorder-dominated phase with a positive droplet exponent ^ > 0, the barrier 
exponent is strictly positive ip > 9 > and leads to some logarithmic dynamics (see Eq |3]) . Since in this paper we 
focus on the value of the barrier exponent ip, we refer the reader to 0, [1, [3] for a detailed description of other 
properties of the droplet dynamics. 

C. Arguments in favor of the value ^ = ds/2 

From the point of view of the dynamical exponent tp, there has been a long-standing difference between 
(i) the directed polymer model, where the assumption that the barrier exponent coincides with the droplet exponent 
has been made from the very first article that has introduced the model [l^ (see section IIIII for more details) 

Usual Assumption for the directed polymer : ip = 9 (14) 

(ii) spin-glasses, where it has been quickly clear that the barrier exponent in strictly bigger than the droplet exponent 
(see section Ivl for more details), because they are distinct below the lower critical dimension. In dimension d = 1, the 
exact solution [1] yields 

ID Spin - glass : -0 = > 6* = -1 (15) 
and in dimension d = 2 these two exponents do not have the same sign 

2D Spin - glass : ^ > > 61 (16) 

The usual explanation of this difference between the two models (i) and (ii) is that the directed polymer case would 
be much more 'simple' than the spin-glass case, that its phase space would be characterized by a single exponent 9, 
whereas in spin-glasses the barrier scaling is not related to the scaling of the free-energy minima. In this paper, we 
propose another scenario, based on the observation that in any disorder system presenting a positive chaos exponent 

= ds/2 — 9 > (and in particular for the directed polymer), the description of the phase space requires at least two 
exponents which are the droplet exponent 9 for free-energy fluctuations (Eq[H]) and the exponent ds/2 that governs 
energy and entropy fluctuations fEg llOp. The assumption of Eq. [T3]is then equivalent to the very strong requirement 
that the global free-energy optimization of order that results from a near cancellation of much bigger energy and 
entropy random contributions of order L'^*/^ is satisfied all along the dynamical trajectories . The alternative scenario 
that we propose in this paper is that for any dynamics containing only local moves of the polymer or interface, the 
barrier exponent ip is equal to the energy-entropy fluctuation exponent ds/2 fEq llOp 

^ = y (17) 

The physical interpretation is that the equality ip = 9 would be possible only via a non-local dynamics that would 
allow a global reorganization of the polymer or interface at each time step, whereas any local dynamics will see barriers 
that are dominated by small-scale contributions of random sign. 

In the remaining sections, we discuss this conjecture for various disordered models and compare with the available 
numerical results. 



III. NON-EQUILIBRIUM DYNAMICS OF DIRECTED POLYMERS IN RANDOM MEDIA 

A. Reminder on the statics 

The directed polymer in a random medium (see for a review) is a model where the various statements of the 
droplet scaling theory have been successfully tested. The exponent 9 of Eq. [6] is exactly known in one-dimension 
9{d = 1) = 1/3 fldi [itI. [TsII and for the mean- field version on the Cayley tree 9{d = oo) = [l^. In finite dimensions 
d = 2,3,4,5,..., the exponent 9{d) has been numerically measured, with values of order 9{d = 2) = 0.244 and 
9{d = 3) = 0.186 [10, l2ll |. The statistics of Eq. [5] for the low-energy excitations as a function of their size / describes 
very well the numerical data in the regime 1 ^ Z <C i in dimensions d = 1, 2, 3 [l^]. Finally, the scaling of Eq. [TOlfor 
the energy and entropy fluctuations have been numerically checked in various dimensions in [i,!!!, 24]. Let us stress 
again that the difference between free-energy fluctuations of Eq. [5] and energy fluctuations of Eq. [TO] can be seen only 
for sufficiently large scale L [1, [H, [l^] (see the discussion before Eq [T^ . 
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B. Discussion of the conjecture ijj = ds/2 = 1/2 for the barrier exponent 

To the best of our knowledge, all papers discussing the non-equilibrium dynamics of the directed polymer seem to 
have assumed the equality of Eq. [T3] between the barrier exponent -0 and the droplet exponent 9. This assumption 
was first made in the very first paper [T^ introducing the directed polymer model, i.e. before the drop let analysis of 
the model Q. More recently, many papers consider that the equality of Eq. [14] has been 'proven' in |25l | up to possible 
logarithmic corrections. In our opinion, the arguments contained in j25l | are plagued by the fact that the authors of 
[25| seem to be unaware of the crucial difference in scaling between free-energy fluctuations and energy fluctuations 
fEqsl^l [TU)) . For instance, they state that the dynamics is controlled by 'energy barriers', which have the 'same scaling 
as free-energy fluctuations', because it is a 'zero-temperature fixed point', but as recalled above, within the droplet 
theory, the properties of the 'zero-temperature fixed point' are instead the different scalings of Eqs [9l and [TOl 

We have explained above in section Hi CI our general arguments in favor of the value -0 = ds/2. For the directed 
polymer of dimension = 1 in a random medium of dimension 1 + d, this corresponds to 

iJDP='^ = \ (18) 
We now compare with the available numerical results. 



C. Comparison with available numerical results on the non-equilibrium dynamics 

Although aging effects for the directed polymer have been fitted with algebraic time scalings by various authors 
[2^ , the more recent work of Kolton, Rosso and Giamarchi [27] shows that 

(a) the growing length L{t) cannot be fitted by a power-law L{t) '-^ t^^^ at large times, although the short time 
relaxation could be fitted with some effective exponent z{T) that strongly depends on temperature. This could explain 
why the first numerical fits [26J see apparent algebraic aging forms. 

(b) the growing length L(t) can be fitted with the logarithmic form L{t) ~ (Int)^/''' at large times, and the value 
of the barrier exponent ij) is asymptotically size and time independent as it should. 

(c) the value of measured in p7j is V' — 0.49 in the last three decades. The interpretation of the authors of [l^] 
that believe in the identity ijj = = 1/3, is that barriers contain strong logarithmic corrections B{L) ^ i^/'^(lni)^. 
Our interpretation is on the contrary that the measured value ~ 1/2 is actually the correct one. 



IV. NON-EQUILIBRIUM DYNAMICS IN DISORDERED FERROMAGNETS 

A. Numerical results on coarsening in disordered ferromagnets 

The non-equilibrium dynamics in pure ferromagnets in the low-temperature phase T < Tc is well understood via the 
characterization of domain coarsening [28| . In the presence of quenched disorder however, the large time behavior of 
the characteristic length scale R(t) of the coarsening process has remained controversial between logarithmic behavior 

R{t) - {\iitf (19) 

with a universal exponent x discussed below, and power-law growth [3l| 

R{t) - (20) 

with an exponent z{T^ e) that depends both on the temperature T and on the disorder strength e. 

However, as stressed in the recent work j33 [. the numerical simulations of coarsening in disordered ferromagnets do 
not really reach the large-scale large-time regime, since the maximal size Rmax measured is sometimes only of order 
Rmax ~ 7 in unit of lattice spacings at the end of the simulation (see for instance Fig 4a of pi] or Fig. 1 of 32]). As 
a consequence, the available numerical simulations on coarsening with quenched disorder are not very conclusive for 
the asymptotic regime of R{t). 
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B. Relation between the exponent x in d — 2 and the directed polymer barrier exponent tpDP 

Since the directed polymer model discussed in prev ious section has been precisely introduced as a model of domain 
wall in two-dimensional disordered ferromagnet [14| . one expects some relation between the barrier exponent ipop 
of the directed polymer in 1 + 1 and the exponent x of Eq 1191 go verning the domain growth in the two-dimensional 
disordered ferromagnet. The first possibility would be simply [30|. 

'^simple — ~i (^-^) 

Vdp 

meaning that the dynamics is governed by the barriers associated to the domain scale R{t). However, Huse-Henley 
have proposed another scenario leading to the higher value 

(2-C) 

XHH = — (22) 

Vdp 

where ^ = 2/3 is the roughness exponent of the directed polymer in 1-1-1. The argument leading to the value of Eq 
can be summarized as follows [3, H^, H^l . The relevant interfaces during the coarsening process are not directed 
polymers but curved polymers with a typical curvature radius of order R{t) itself. It can be consider as directed up 
to the size / where the roughness is of the same order of the curvature P/R(t) yielding l{t) ^ {R{t)y^^^~^\ The 
barriers associated to these directed parts scale as (^(i))'^"^ - leading to Eq. [H 

Note that these arguments usually go along with the assumption ipop = 0dp — 1/3 (see previous section on 
the directed polymer) yielding the values Xsimpie = 3 and xhh = 4 [13, [H, H^]- With the value ipDP — 1/2 of 
our conjecture discussed in previous section for the directed polymer, the values of the exponent x are respectively 
Xsimpie = 2 and Xhh — 8/3. In the following, we argue that within our analysis, it is the value x simple = 2 which is 
natural for disordered ferromagnets in d = 2. 

C. Conjecture ijj — ds/2 — {d~ l)/2 for the barrier exponent 

Within our analysis where the barrier exponent -0 is governed by the dimensionality ds of the interface that 
determines the energy and entropy fluctuations (Eq llOp . we expect that in disordered ferromagnets where ds = d — 1, 
the barrier exponent is 

^ = y-^ (23) 

irrespectively of the directed or curved nature of the interface, since it is governed by small-scale contributions. In 
particular, in domain coarsening, we expect the 'simple' relation that generalizes Eq. [2T] 

Again, as explained above, numerical data on coarsening with disorder do not allow a precise measure of the exponent 
X because R{t) of Eq. [12] is never very large in simulations (see [s^l for a more detailed discussions of the results of 
various fits). 

V. NON-EQUILIBRIUM DYNAMICS IN SPIN-GLASSES 

Many numerical works have studied non-equilibrium properties in spin-glasses. Here again, there is a controversy 
between logarithmic dynamics (see for instance [H, [s^]) and algebraic dynamics (see for instance [ssl. [36|). 

A. Discussion of the conjecture = ds/2 for spin-glasses in d = 2 

In dimension d — 2, there is no spin-glass phase because the droplet exponent 9 is negative 6* < so that Tc — 0. 
Nevertheless, it is interesting to measure the values of the droplet exponent 9 and of the fractal dimension ds of the 
surface of excitations above the ground state. Recent estimates are 9 ~ —0.287(4) (see [13 and references therein) 
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and ds — 1.274(2) (see [ssl] and references therein). Note that it has been recently argued that these interfaces are 
described by SLE evolutions implying some simple relation between 9 and dg [soj . 

Using ds — 1.274 [sH], the present conjecture ijj — ds/2 for the barrier exponent would corresponds to a numerical 
value of order 

2D Spin - glass : i; = ^ - 0-637 (25) 

In a recent work, Amoruso, Hartmann and Moore [40| have tried to measure the barrier exponent of the highest 
barrier of systems of sizes L < 40 yielding the numerical bounds 0.25 < ip < 0.54. The uncertainty shows the difficulty 
of the numerical measure of so that their upper bound does not seem to us sharp enough to rule out the value of 

B. Discussion of the conjecture = ds/2 for spin-glasses in d = 3 

In dimension d = 3, the droplet exponent 9 is positive 9 > so there exists a spin- glass phase with Tc > 0. Recent 
estimates for the droplet exponent 9 and for the fractal dimension ds of the surface of excitations above the ground 
state are respectively 9 ~ 0.19(2) (see [4l| and references therein) and ds — 2.6 (see and references therein). Using 
the latter, the present conjecture ip = ds/2 for the barrier exponent would correspond to a numerical value of order 

3D Spin - glass : V = y =i 1-3 (26) 

This has to be compared with the value i/j ~ 1.0 estimated by Berthier and Bouchaud from their aging simulations 
[33 |. Again, the precision of this numerical estimate does not seem sufficient to rule out the value of Eq[26l We refer 
the reader to [3J] for the experimental values of the exponent ip reported in the literature, that varies between 0.3 
and 1.9 m. 

VI. CONCLUSION 

In this paper, we have proposed that, in disordered systems characterized by a positive chaos exponent ( = ds/2— 9 > 
0, the large time dynamics is governed by the barrier exponent ip = ds/2. We have explained why this value ip = ds/2 
is natural within the droplet scaling picture, where the exponent ds/2 governs the energy-entropy fluctuations (Eq [TO|) 
and is greater than the droplet exponent 9 of free-energy fluctuations (Eq[21). We have then discussed our conjecture 
for the following models 

(i) for the directed polymer where ds = 1, our conjecture gives ip = 1/2 in all dimensions ; 

(ii) for disordered ferromagnets where ds — d — 1, our conjecture corresponds to ip = {d — l)/2 

(iii) for spin-glasses where interfaces have a non-trivial fractal dimension ds known numerically, our conjecture ip — 
ds/2 gives numerical predictions in d = 2 and d = 3. 

In each case, we have compared with the available numerical data on -0, in particular with the work of Kolton-Rosso- 
Giamarchi[23| who have measured the barrier exponent ip ~ 0.49 for the non-equilibrium dynamics of a directed elastic 
string. For disorder spin models, either disordered ferromagnets or spin-glasses, the available numerical estimates of 
Ip are not sufficiently precise to support or exclude our conjecture. 

If the conjecture ip — ds/2 is correct, this means that the numerical measure of this barrier exponent in dynamical 
simulations requires to study samples of sizes L sufficiently large, where in the corresponding statics, the free-energy 
fluctuations and the energy-entropy fluctuations have reached their asymptotic regimes of Eqs [TUl and [91 i.e. one needs 
to be in the regime of Eq[T2]for the statics. We hope that this explicit static criterion will help to identify the regime 
where dynamic simulations are likely to measure the asymptotic barrier exponent ip relevant at large scales. 
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